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ABSTRACT 


The  Cagniard-de  Hoop  method  is  ideally  suited  to  the  analysis  of  wave  propagation 
problems  in  stratified  media.  The  method  applies  to  the  integral  transform  represen¬ 
tation  of  the  solution  in  the  transform  variables  (s.p)  dual  of  the  time  and  transverse 
distance.  The  objective  of  the  method  is  to  make  the  p-integral  take  the  form  of  a 
forward  Laplace  transform,  so  that  the  cascade  of  the  two  integrals  can  be  identified 
as  a  forward  and  inverse  transform,  thereby  making  the  actual  integration  unnecessary. 
That  is,  the  exponent,  —sw(p)  is  set  equal  to  — sr,  with  r  varying  from  some  (real) 
finite  time  to  infinity.  As  usually  presented,  the  p-intcgral  is  deformed  onto  a  contour 
on  which  the  exponent  is  real  and  decreases  to  — oe  as  p  goes  to  infinity.  We  have  found 
that  it  is  often  easier  to  introduce  a  complex  variable  r  for  the  exponent  and  carry  out 
the  deformation  of  contour  in  the  complex  r-domain.  In  the  r-domain  the  deformation 
amounts  to  “closing  down”  the  contour  of  integration  around  the  real  axis  while  taking 
due  account  of  singularities  off  this  axis. 

Typically,  the  method  is  applied  to  an  integral  that  represents  one  body  wave  plus 
other  types  of  waves.  In  this  approach,  the  saddle  point  of  w(p)  that  produces  the  body 
wave  plays  a  crucial  role  because  it  is  always  a  branch  point  of  the  integrand  in  the 
r-domain  integral.  Furthermore,  the  paths  of  steepest  ascent  from  the  saddle  point 
are  always  the  tails  of  the  Cagniard  path  along  which  w(p)  —>  oc.  That  is,  the  image 
of  the  pair  of  steepest  ascent  paths  in  the  p-domain  is  a  double  covering  of  a  segment 
of  the  Re  r  axis  in  the  r-domain.  The  deformed  contour  in  the  p-domain  will  be  only 
the  pair  of  steepest  ascent  paths  unless  the  original  integrand  had  other  singularities 
in  the  p-domain  between  the  imaginary  axis  and  this  pair  of  contours.  This  motivates 
the  definition  of  a  primary  p-domain — the  domain  between  the  imaginary  axis  and  the 
steepest  descent  paths — and  its  image  in  the  r-domain — the  primary  r-domain.  In 
terms  of  these  regions,  singularities  in  the  primary  p-domain  have  images  in  the  pri¬ 
mary  r-domain  and  the  deformation  of  contour  onto  the  real  axis  in  the  r-domain  must 
include  contributions  from  these  singularities. 

This  approach  to  the  Cagniard-de  Hoop  method  represents  a  return  from  de  Hoop’s 
modification  to  Cagniard's  original  method,  but  with  simplifications  that  make  the 
original  method  more  tractable  and  straightforward.  This  approach  is  also  reminiscent 
of  van  der  Waerden’s  approach  to  the  method  of  steepest  descents,  which  starts  exactly 
the  same  way.  Indeed,  after  the  deformation  of  contour  in  the  r-domain,  the  user  has 
the  choice  of  applying  asymptotic  analysis  to  the  resulting  “loop”  integrals  (Watson’s 
lemma)  or  continuing  to  obtain  the  exact,  although  usually  implicit,  time  domain  so¬ 
lution  by  completing  the  Cagniard-de  Hoop  analysis. 

In  developing  the  method  we  examine  the  transformation  from  a  frequency  domain 
representation  of  the  solution  (w)  to  a  Laplace  representation  (s).  Many  users  start 
from  the  frequency  domain  representation  of  solutions  of  wave  propagation  problems. 
There  are  issues  of  movement  of  singularities  under  the  transformation  from  u  to  s  to 
be  concerned  with  here.  We  discuss  this  anomaly  in  the  context  of  the  Sommerfeld 
half-plane  problem. 
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INTRODUCTION 

The  Cagniard-de  Hoop  method,  or  the  Cagniard  (1939,  19G2)  method  as  modified 
by  de  Hoop  (1958),  is  ideally  suited  to  the  analysis  of  wave  propagation  problems  in 
stratified  media.  The  method  applies  to  the  integral  transform  representation  of  the 
solution  in  the  transform  variables  (s,p)  dual  of  the  time  and  transverse  distance, 
typically  of  the  form, 

u{r,z,s)  =  J^f(p)exp{-sw{p,r,z)}dp ,  (0.1) 

with  w  having  the  form, 

j  _ 

w(p,  r,z)  =  pr  +  Y,  hj  sjp]  -  p2-  (0.2) 

j= i 

In  this  equation  each  hj  represents  the  total  vertical  travel  distance  in  layer  j,  includ¬ 
ing  multiple  reflection  path  lengths  on  the  trajectory  from  the  source  to  an  observation 
point,  (a;,  ^).  In  a  two-dimensional  problem  ( x ,  z),  r  —  |x|  denotes  the  transverse  range 
and  pj  —  l/cj  is  the  slowness  in  that  layer.  In  a  three-dimensional  problem  (x,y,z), 
r  =  \/x^  +  y-  and  p  is  the  component  of  the  slowness  vector  which  is  colinear  with  the 
vector,  x  =  ( x ,  y ).  Furthermore,  p2  =  l/c2  +  q 2,  with  q  the  component  of  the  slowness 
vector  in  the  orthogonal  direction  to  x.  The  contour  of  integration  T  is  the  imaginary 
p-axis.  The  integral  (0.1)  is  also  known  as  a  generalized-ray  wave  constitutent  or, 
more  simply,  as  a  generalized  ray  (Spencer,  1960;  Cisternas,  et  al.,  1973). 

The  objective  of  the  method  is  to  make  the  slowness  p-integral  take  the  form  of 
a  forward  Laplace  transform  in  time  (possibly,  a  sum  of  such  transforms).  Each  of 
these  will  be  of  the  form, 

/oo  dry 

/(p(r))^T  exp{— sr}  dr.  (0.3) 

If  this  is  accomplished,  then  the  cascade  of  this  integral  with  the  inverse  Laplace 
transform  from  .s  back  to  time  t  can  be  identified  as  a  forward  and  inverse  transform 
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pair,  making  the  actual  inversion  integration  unnecessary.  To  achieve  this,  the  contour 
T — the  imaginary  axis-  -is  deformed  onto  a  path  of  integration  on  which  w{p)  in  (0.1) 
is  real  and  approaches  infinity.  Then  w(p)  is  set  equal  to  a  real  variable  r  that  varies 
from  some  finite  time  to  infinity. 

Excellent  expositions  of  the  Cagniard-de  Hoop  method  may  be  found  in  the  exten¬ 
sive  papers  of  A.  T.  De  Hoop,  including  his  thesis  (1958)  and  his  recent  tutorial  (1988); 
also,  Aki  and  Richards  (1980)  and  van  der  Hijden  (1987)  are  important  sources. 

The  presentation  of  the  method  often  begins  with  a  “folding  over”  of  the  contour 
so  that  one  endpoint  of  the  integral  is  at  p  =  0.  This  leads  to  a  representation 
for  u(r,  z,s)  as  the  imaginary  part  of  an  integral  in  r.  In  this  approach  the  p-path 
contains  a  section  of  contour  starting  from  p  =  0  where  the  integrand  is  real,  yielding 
a  zero  contribution  to  the  final  result.  However,  this  one-sided  representation  is  not 
essential  to  the  method  and  we  will  work  with  the  original  two-sided  representation 
allowing  the  upper  and  lower  contributions  to  combine  at  the  end  of  the  analysis,  in 
particular  canceling  away  the  initial  segment  just  discussed.  Van  der  Hijden  is  an 
exception;  he  also  applies  the  Cagniard-de  Hoop  method  to  the  integral  with  both 
endpoints  at  infinity  as  we  do  here,  as  does  de  Hoop  on  occasion. 

We  have  found  that  allowing  r  to  have  equal  status  with  p  as  a  complex  variable 
and  considering  both  these  complex  planes,  as  Cagniard  did  [and  also  Dix  (1954) 
and  Garvin  (1956),  following  Cagniard],  is  an  important  analytic  tool.  Thus,  we  take 
full  advantage  of  the  standard  p-domain  analysis  in  locating  the  singularities  of  the 
integrand  in  (0.1),  but  we  carry  out  the  deformation  of  contour  in  the  r-domain, 
where  it  amounts  to  “closing  down”  the  contour  of  integration  around  the  real  axis 
while  taking  due  account  of  singularities  off  this  axis.  This  is  often  far  easier  than 
finding  the  curves  where  only  R ew  varies,  or  equivalently,  the  curves  where  Ini  to  is 
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constant.  This  deformation  of  contour  is  further  facilitated  by  having  both  endpoints 
at  infinity,  rather  than  having  one  (unphysical)  endpoint  at  the  origin  in  p. 

In  this  approach  the  saddle  points  of  the  exponent  in  the  original  integral  play 
a  crucial  role  because  they  are  always  branch  points  of  the  integrand  f  {p{r))dp  /  dr , 
in  the  r-domain  integral.  On  the  other  hand,  we  will  show  that  the  images  of  the 
branch  points  at  ±pj  in  equation  (0.2)  are  points  of  analyticity  of  dp/dr  in  the  r- 
domain.  Furthermore,  branch  points  of  the  amplitude  f(p)  of  the  same  order  as 
in  the  exponent  at  any  pj  are  regular  points  of  the  amplitude  in  (0.3)  in  the  r- 
domain.  Other  singularities  of  /(p(r))  will  provide  corresponding  singularities  in  the 
r-domain.  Thus,  it  is  the  saddle  points  and  these  other  singularities  that  provide 
the  important  structure  of  the  integrand  in  the  r-domain.  The  exact  evaluation  of 
the  integral  involves  residues  from  the  poles  enclosed  in  the  deformation  process  plus 
loop  integrals  around  each  of  the  branch  points  in  the  r-domain  that  are  encountered 
in  the  deformation  process.  It  is  the  latter  type  of  integral  that  leads  to  integrals  of 
the  form  (0.3);  the  inverse  transform  of  residues  at  the  poles  is  more  straightforward. 
The  deformation  of  contour  requires  the  estimate  of  the  integral  on  a  circular  arc 
“at  infinity.”  The  estimate  that  this  integral  approaches  zero  with  increasing  radius 
follows  from  Jordan’s  lemma. 

A  particular  saddle  point,  p,,  located  on  the  Rep-axis  plays  an  especially  impor¬ 
tant  role.  Indeed,  the  essence  of  our  point  of  view  can  be  explained  in  the  context 
of  this  special  saddle  point  with  the  aid  of  Figures  0.1  and  0.2.  In  the  former,  the 
contour  S  is  made  up  of  two  steepest  ascent  1  paths  from  the  saddle  point  ps  when 
r  ^  0.  We  also  show  the  contour  T.  We  will  refer  to  the  region  between  T  and  S  as 
the  primary  p- domain.  It  will  be  shown  that  the  image  T'  of  T  is  the  hyperbolic-like 
curve  of  Figure  0.2  and  the  image  S'  of  the  contour  S  is  a  semi-infinite  segment  of 

'Paths  of  steepest  ascent  and  descent  are  paths  along  which  Retc  changes  maximally;  they  are 
also  paths  on  which  Im  ?r  is  constant. 
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Fig.  0.1.  Path  T  and  steepest  descent  path  for  an  example  with  J  —  2,  h\  =  h2  =  1.2, 
r  =  2.5,  1/pi  =  3,  l/p2  =  4. 


Fig.  0.2.  Image  contours  for  example  of  previous  figure. 
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the  real  axis,  doubly  covered,  starting  from  ra,  the  image  of  the  saddle  point  p3.  This 
implies  that  the  two  steepest  ascent  paths  which  comprise  S  are  each  at  least  a  tail 
of  a  Cagniard  path.  For  r  0,  ts  is  to  the  right  of  the  point  where  T'  crosses  the 
Rer  axis.  The  lower  part  of  S  maps  onto  the  real  axis,  “below”  a  branch  cut  from  ts 
to  oc  with  arg(r«  —  r)  =  ix  here;  the  upper  part  of  S  maps  onto  the  real  axis,  “above” 
a  branch  cut  from  ts  to  oo  with  arg(r5  —  r)  =  —  n  here.  The  image  of  the  primary 
p-domain  is  the  region  between  T'  and  S'  in  the  r  domain.  We  will  call  this  region 
the  primary  r-domain.  We  remark  that  in  this  region  |  arg(r,  —  r)|  <  tt. 

Finding  the  p-path  along  which  w(r,  z,p)  is  real  and  approaches  +oo,  is  equivalent 
to  closing  the  contour  T'  down  around  the  real  axis  in  the  r-domain.  If  there  are  no 
complex  singularities  of  the  integrand,  f  (p(r))dp  /  dr ,  in  the  primary  r-domain  then 
we  replace  F'  by  S'.  The  integrals  above  and  below  the  cut  are  two  integrals  of  the 
tvpe  (0.3),  with  two  different  branches  of  the  function  p(r)  taken  in  each  integral. 
Typically  this  amounts  to  taking  one  or  the  other  sign  of  a  square  root.  In  this 
case,  the  two  steepest  ascent  paths  comprising  S  are  the  Cagniard  paths  for  two 
semi-infinite  integrals  whose  form  in  the  r-domain  is  exactly  (0.3).  From  the  method 
used  to  arrive  at  the  deformed  paths  as  segments  of  S',  it  is  clear  how  to  deal  with 
the  singularities  on  the  real  r  axis:  the  contour  of  integration  passes  above  such  a 
singularity  on  the  upper  part  of  the  path  S'  and  passes  below  such  a  singularity  on 
the  lower  part  of  the  path  S'. 

If  there  are  singularities  in  the  primary  r-domain,  then  one  must  take  account 
of  them  in  carrying  out  the  deformation  of  contour.  However,  singularities  in  the 
primary  r-domain  are  the  images  of  singularities  in  the  primary  p-domain.  Thus, 
we  can  find  these  singularities  in  the  p- plane,  which  is  often  easier.  Furthermore, 
we  can  be  sure  that  singularities  to  the  left  of  the  primary  p-domain  are  also  to 
the  left  of  the  primary  r-domain  and  have  no  effect  on  the  contour  deformation 
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process.  Somewhat  similarly,  singularities  to  the  right  of  the  primary  p-domain  are 
also  outside  the  primary  r-domain.  They  lie  in  the  region  where  |arg(r,  -  r)|  >  7r, 
usually  described  as  being  on  a  “second  Riemann  sheet”  of  the  r-domain.  While  they 
do  not  obstruct  the  deformation  of  contour  from  T'  to  S',  they  can  be  close  enough 
to  S'  that  their  “influence”  is  felt  in  the  integration  along  S'.  Examples  of  this  type 
are  discussed  in  the  above  cited  literature. 

We  have  twice  included  the  caveat  r  /  0.  The  case  r  =  0  is  degenerate;  T  =  S  and 
the  primary  p-domain  collapses  to  the  imaginary  axis  with  its  image  being  the  right 
half  Rer  axis,  starting  from  t3  =  r(p  =  0),  doubly  covered  just  as  in  Figure  0.2.  Here, 
there  is  no  chance  for  other  singularities  to  be  “caught”  between  the  tw’o  contours 
and  the  entire  field  arises  as  the  body  wave  contribution  due  to  the  saddle  point. 
(Note  that  this  is  the  case  of  a  vertical — possible  multiply  reflected — ray  trajectory.) 
As  above,  however,  singularities  to  the  right  of  T  =  S  may  yet  appear  close  enough 
to  S'  to  influence  the  value  of  the  integral  over  S'. 

In  summary,  our  approach  may  be  viewed  as  an  attempt  to  replace  the  contour 
T'  by  the  contour  S'.  When  there  are  no  singularities  in  the  primary  r-domain- 
equivalently,  no  singularities  in  the  primary  p-domain-the  replacement  is  justified  by 
Cauchy’s  theorem;  when  there  are  singularities  in  the  primarj'  r-domain,  we  must 
account  for  these  singularities  in  the  contour  deformation  and  S'  is  only  the  tail  of 
(one  of  the)  loop  contours  in  the  r-domain. 

This  approach  to  the  Cagniard-de  Hoop  method  is  reminiscent  of  van  der  Waer- 
den’s  (1951)  approach  to  the  method  of  steepest  descents,  which  starts  exactly  the 
same  way.  Indeed,  after  the  deformation  of  contour  in  the  r-domain,  the  user  has  the 
choice  of  (i)  applying  asymptotic  analysis  to  the  resulting  “loop”  integral(s)  [Wat¬ 
son's  lemma]  (Bleistein,  1984)  or  (ii)  continuing  to  obtain  the  exact,  although  usually 
implicit,  time  domain  solution  by  fomoleting  the  Cagniard-de  Hoop  analysis.  In 
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the  former  case,  one  obtains  the  large  s  expansion,  from  which  a  progressing  wave 
expansion  in  the  time  domain  can  be  derived. 

In  some  sense,  we  have  come  full  circle  back  to  Cagniard’s  (1939,  1962)  approach, 
using  two  complex  planes.  In  that  work  the  author  described  the  mapping  from  p 
to  r  in  much  more  detail  than  we  think  is  necessary.  Indeed,  de  Hoop  (1958)  points 
out  that  one  of  his  motivations  for  his  modification  of  the  Cagniard  method  was  the 
difficulties  of  coping  with  the  complex  mapping.  Furthermore,  the  important  role 
that  we  see  for  the  path  of  steepest  ascent  was  not  a  part  of  that  original  method. 
We  believe  that  our  “equal  status”  point  of  view  extracts  the  best  of  both  methods. 


In  the  next  section,  we  discuss  the  transformation  from  Fourier  transforms  in 
traditional  wave  number  k  (one  transverse  dimension)  or  wave  vector  k  and  frequency 
u j  to  slowness  p  or  slowness  vector  p  and  Laplace  transform  vaiiable  s.  We  do  this  for 
two  reasons.  First,  for  some  users  deriving  a  solution  representation  in  the  former  pair 
of  variables  is  more  natural.  Second,  in  the  former  variables  we  can  treat  a  traditional 
noncau'al  probL...  the  Soramcifdd  problem  with  acute  angle  of  incidence  of  the 


plane  wave — that  is  not  amenable  to  analysis  in  the  Laplace  domain.  In  the  third 
section,  we  discuss  the  transformation  from  p  to  r  with  emphasis  on  features  that  we 
believe  are  important  to  our  approach  to  ihe  method.  Finally,  we  appl}  our  approach 
to  some  standard  examples  of  the  Cagniard  method  and  then  to  the  nonstandard 
Sommcrfeld  example. 


TRANSFORMATION  OF  FOURIER  INTEGRALS 

In  our  analysis  of  waves  in  layered  media  we  start  from  a  Fourier  transform  in  time, 
rather  than  from  tin;  Laplace  transform  typical  of  the  literature  of  the  Cagniard-de 
Hoop  method.  We  think  this  is  a  useful  choice  because  Fourier  representations  for 
solutions  of  wave  propagation  problems  are  far  more  common  than  Lapla^  icpresen- 
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tat  ions  in  the  modern  literature  with  the  Cagnaird-de  Hoop  method  being  almost  the 
sole  exception  to  this  trend. 

For  "causal  "  problems,  the  Fourier  representation,  just  as  the  Laplace  representa¬ 
tion,  is  a  transform  on  a  semi-infinite  domain  from  some  finite  time  to  infinity.  Most 
problems  of  interest  start  from  some  finite  time.  Indeed,  even  some  problems  de¬ 
scribed  in  terms  of  a  plane  wave  incident  from  infinity  can  be  reformulated  as  causal 
problems.  We  will  discuss  this  point  further  below  in  connection  with  the  Sommerfeld 
half-plane  problem.  Thus,  for  a  function  V(x,z,t),  we  define  in  the  causal  case 

v(x,z,u)  =  /  V  (.t,  z ,t)  exp l  iuit }  dt .  (0.4) 

Jin 

It  follows  that  the  temporal  Fourier  transform  is  an  analytic  function  of  uj  in  some 
upper  half  m-plane,  usually  not  including  the  real  axis  where  singularities  of  the 
transform  reside.  Of  course,  this  upper  half-plane  is  exactly  the  right  half-plane  of 
analyticity  ot  mo  Laplace  transform.  The  value  of  the  solution  representations  for 
real  values  of  uj  are  oVaii  y[  analytic  continuation  from  this  upper  half-plane. 
In  fact,  we  can  and  v. ih  consider  only  the  upper  right  half-plane  Rem  >  0,  and 
determine  the  solution  elsewhere  by  analytic  continuation.  With  this  convention,  the 
specific  branches  of  multi-valued  functions  of  m  and  the  spatial  transform  variable(s) 
are  uniquely  defined,  as  is  the  avoidance  of  these  singularities  on  the  contour(s)  of 
integration  of  the  inverse  transform. 

In  most  cases  of  interest  in  wave  propagation  problems,  the  lower  boundary  of  the 
half-plane  of  analyticity  is  the  Re  m  axis.  We  assume  that  is  the  case  here,  although 
it  is  not  crucial  to  the  analysis  below.  However,  it  does  allow  us  to  define  the  quarter 
plane  of  our  attention  by  the  condition 

w  --  |mj  exp  {’O'}  0  <  a  =  argm  <  7r/2.  (0.5) 
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For  problems  in  two  spatial  dimensions,  the  Fourier  (A\u;)  representation  to  be 
studied  takes  the  form 

r(j\  =  /  g(k,^;)cxp{ip(k\^\.i\z)}  dk\ 

J  —  O C 

(0.6) 

j  _ 

r,z)  =  kT  +  ^hjy/^/cj)2  -  k'2. 

j= i 

The  function  /;  is  assumed  to  have  the  structure 

g(k^)  =  (^)ng(k/^,  1).  (0.7) 

This  is  the  form  of  g  for  a  point  source:  for  more  complex  sources,  the  solution  is 
obtained  by  convolution  with  the  solutions  of  the  type  presented  here. 

With  our  assumption  that  Inu'  >  0.  there  are  no  singularities  of  p  (or  g)  on  the 
Re  /.-axis,  the  path  of  integration  in  (0.6).  In  particular,  the  branch  points  at  uJ 
lie  in  the  first  quadrant  (or  positive  imaginary  axis)  of  the  A-plane,  and  the  branch 
points  at  lie  in  the  third  quadrant  (or  negative  imaginary  axis),  because  we 

have  fixed  Re_c  >  0. 

Let  us  introduce  the  change*  of  variable  of  integration 

k  —  p ,  (0.8) 

and  consider  the  image  of  the  integration  contour  in  the  p-plane  and  the  location  of 
singularities  of  the  integrand  there,  as  well. 

First,  note  that  for  k  real  and  argu,*  =  o  (0.5).  on  the  path  of  integration 

argp=—  a.  (0.9) 

This  means  that  as  arg  v  increases  from  zero  to  ir/2  the  contours  of  integration  are 
a  family  of  straight  lines  through  the  origin,  passing  from  the  first  quadrant  to  the 
third. 
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To  achieve  the  form  (0.1),  it  is  necessary  to  pick  a  =  tt/2.  While  in  principle  one 
could  attempt  an  analog  to  the  Cagniard-de  Hoop  theory  with  other  choices  of  a, 
tins  choice  is  optimal  as  we  will  show  later.  In  this  case,  the  contour  of  integration  is 
tln>  imaginary  axis  starting  at  -Hoc  and  ending  at  —ioc.  Later,  we  will  reverse  this 
direction  by  multiplying  the  integrand  by  —1. 


Before  proceeding  any  further  with  this  line  of  analysis,  let  us  examine  the  singu¬ 
larities  of  tin'  integrand.  In  particular,  note  from  (0.G)  that 

f  ■/  _ 1 


r(A-., 


p.r  +  Y,  h]  \Al/Ti)2 


p- 


(0.10) 


L  ;=1  J 

Thus,  the  integrand  in  p  has  branch  points  at  ±l/c;,  j  —  1, . . . , ./.  and,  as  a  increases 
from  zero  to  z/2,  the  contour  simply  rotates  away  from  these  branch  points. 


On  the  other  hand,  suppose  that  g(k, ce)  had  singularities  in  the  second  or  fourth 
quadrant.  Then  we  would  have  to  take  account  of  the  effect  of  passing  the  contour 
of  integration  "through”  these  singularities.  For  poles,  that  would  simply  amount  to 
including  a  residue  as  part  of  the  transformation  process;  for  branch  points,  a  loop 
integral  enclosing  the  branch  point  would  have  to  be  included  in  the  analysis.  Such 
singularities  arise  only  in  noncausal  problems — an  example  of  this  type  is  included  in 
a  later  section. 


Choosing  a  =  z/2  amounts  to  evaluating  the  Fourier  transform  for  a  purely 
imaginary  value  of  the  transform  variable  u j.  Alternatively,  we  can  set 


u!  —  is 


(0.11) 


with  .s-  real  and  positive  when  o  =  argu;  =  z /2. 

What  we  have  just  done  is  change  a  Fourier  transform  into  a  Laplace  transform. 
However,  we  have  done  it  with  enough  care  to  alert  the  reader  to  the  possibility 
of  contributions  from  singularities  of  the  integrand  of  a  type  we  will  see  below  in 
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a  noncausal  example.  We  note  that  converting  the  forward  Fourier  transform  to  a 
forward  Laplace  transform  implies  that  the  inverse  transform  now  takes  the  form  of 
the  usual  inverse  Laplace  transform. 

We  proceed  below  with  our  analysis  for  s  real  and  positive.  All  results  can  be 
extended  to  complex  s  by  analytic  continuation  techniques;  however,  this  is  not  nec¬ 
essary  for  the  Cagniard-de  Hoop  method  since  the  explicit  inversion  integration  is  not 
performed. 

We  have  one  other  small  trick  to  offer.  We  use  the  transformation  (0.8)  for  x 
positive  and  we  use  k  =  —up  for  x  negative.  That  is,  we  use  k  =  upsgn(x).  This 
reverse.,  the  orientation  of  the  contour,  but  introduces  a  compensating  minus  sign  in 
dk/dp.  These  are  nullifying  effects. 

Furthermore,  functions  of  k2,  such  as  the  square  roots  appearing  in  are  unaf¬ 
fected  by  this  trick,  while  the  expression  kx  is  transformed  into  up\x\.  This  will  prove 
useful  below.  The  amplitude  g  need  not  be  a  function  of  k 2  and  some  care  is  required 
in  dealing  with  this  function. 

In  summary,  taking  account  of  this  additional  trick  and  the  contour  rotation,  as 
well  as  the  change  of  variables  from  u  to  s,  we  set 

k  =  isp  sgn(x).  (0.12) 

For  this  function  dw/dp  >  0  at  p  =  0  for  x  nonzero.  Since  it  is  our  intention  to 
make  t  =  tr  the  new  variable  of  integration,  this  will  insure  that  the  image  T  of  the 
Rep  axis,  the  contour  of  integration  in  (0.1),  has  an  orientation  that  is  independent 
of  sgn(.r).  This  was  the  motivation  for  using  a  change  of  variables  that  depends  on 
sgn(.r). 
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Under  this  combined  change  of  variables  and  rotation  T,  the  contour  of  integration 
in  (0.1)  is  the  imaginary  axis  oriented  from  —ioo  to  zoo  and 

f(p)  =  -9(P  sgn(x),l)  =  ~(is)~ng{ispsgn{x),is), 


u(r,z,s)  —  v(x,  z,is)(is)  "  1 


(0.13) 


Here  we  have  allowed  for  u  to  be  the  true  transformed  field  within  a  power  of  is.  In 
general,  we  neglect  multiplication  by  powers  of  s  here,  since  they  can  be  dealt  with 
back  in  the  space/ time  domain.  Integer  powers  of  s  are  equivalent  to  differentiation 
or  integration  in  the  time  domain;  noninteger  powers  of  s  correspond  to  fractional 
derivatives  or  integrals  in  the  time  domain,  which  are  equivalent  to  convolution  with 
fractional  powers  of  t.  The  extra  minus  sign  in  the  transformation  from  g  to  f  arises 
from  the  fact  that  the  image  of  the  integration  path  in  (0.6)  is  the  imaginary  axis 
oriented  downward ,  whereas  T  in  (0.1)  is  oriented  upward.  Had  we  not  assumed  the 
form  of  g  in  (0.7),  then  we  would  not  have  ended  up  with  /  being  a  function  of  p 
alone  in  this  equation.  Authors  starting  from  Laplace  transforms  assume  this  form 
of  /  to  begin  their  analysis. 


For  three-dimensional  problems,  we  start  from  a  representation  of  the  form 

/OO 

g(k i  ,  A'2,  w)  exp {iip(kuk2,  to,  x ,  y,  c)}  dk{  dk2, 

-OO 

(0.14) 

tp(ki,k2,u,x,y,z)  =  kix  +  k2y  +  J2  -  A-f  -  k\. 

i= i 

Again,  we  assume  that  g  has  the  special  form 


g(k\ ,  k2,uj)  =  {u)n  g{kxjuj,  k2/uj,  1). 


(0.15) 


To  transform  this  integral,  we  first  introduce  the  polar  coordinates  (r,  </>), 

x  —  rcos<t>,  y  —  r  sin</>.  (0.16) 


12 


Bleistein  and  Cohen 


Cagniard  Method,  Complex  Time 


Next,  we  apply  the  de  Hoop  transformation  (de  Hoop,  1960)  in  the  A>domain, 


k  i  =  k  i  cos  d>  —  K2  sin  0, 


k2  =  k i  sin  <p  +  k,2  cos  0; 


(0.17) 


and  we  note  that 


kxx  +  k2y  =  Kir. 


(0.18) 


In  these  new  variables,  (0.14)  becomes 

/oo 

g(Ki,K2,uj)exp{i<p(Ki,K2,u),x,y,z)}  dtci  dn2, 

-OO 


Wi(KUK2,iu,x,y,z)  =  Kir  +  ^2  hjyj(u/cj )2  -  k\  -  k%,  (0.19) 


j= i 


g{Ki,K2,u)  =  g(ki,k2,<jj). 


The  function  g  retains  the  property  (0.15).  We  proceed  with  the  same  scaling  as 
above,  except  that  now  we  must  scale  both  and  k2.  The  rotation  in  u>  rotates 
both  transformed  integrals.  The  same  caveats  above  must  be  observed  here  regarding 
passing  through  singularities. 

In  summary, 

K\  =  isp,  k2  =  sq,  u>  =  is,  (0.20) 


leads  to  the  representation  of  the  form  of  (0.1),  except  that  now  both  /  and  w  depend 
on  q: 

f(p)  =  7,1)  =  ~{is)~ng(isp,sq,is), 


(0.21) 


u(r,z,s)  ~  (is)  "  xs  xv(x,y,  z,is). 


13 


Bleistein  and  Cohen  Cagniard  Method,  Complex  Time 

For  the  purposes  of  the  Cagniard-de  Hoop  method,  we  treat  q  as  a  parameter  whose 
presence  need  not  be  specifically  noted  in  the  arguments  of  these  functions.  However, 
it  should  be  noted  that  q  changes  the  meaning  of  the  slowness  value  l/pj,  in  that 


Pj  =  1  /cj,  two  dimensions, 


Pi  = 


three  dimensions. 


(0.22) 


Furthermore,  after  applying  the  Cagniard  method  to  determine  U(r,  z,t),  it  is  neces¬ 
sary  to  integrate  the  result  with  respect  to  q  over  (  —  00,00). 

In  summary,  we  have  shown  how  Fourier  representations  of  both  two-  and  three- 
dimensional  wavefields  in  stratified  media  can  be  reduced  to  the  same  integral  form, 
Equation  (0.1). 

ANALYSIS  OF  THE  EXPONENT 

Here  we  consider  the  complex  change  of  variables 

r  =  w{p,r,z)  =  pr  +  Y,hj\/P2i  ~  P2,  (0.23) 

i= 1 

from  p  to  t.  Under  this  change  of  variables,  the  integral  in  (0.1)  becomes 

u(x,  z,s )  =  f{p{r))^e~ST  dr.  (0.24) 


We  need  to  know  the  image  T'  in  the  r-plane  of  the  contour  of  integration  in  (0.1). 
That  original  contour  T  is  just  the  imaginary  axis.  Because  we  want  to  deform  this 
image  contour  onto  the  Re  r  axis,  we  also  need  to  understand  the  conformal  mapping 
of  the  image  region  between  T'  and  the  Re  r  axis.  In  this  manner,  we  will  identify 
those  singularities  that  must  be  crossed  in  deforming  one  contour  into  the  other. 
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The  Image  of  f 


We  consider  now  the  image  of  the  contour  T  under  the  mapping  (0.23).  For  this 
image,  let  us  set 

p  =  io ,  —  oc  <  o  <  oo.  (0.25) 


With  this  choice, 


Re  r  =  ^  hj  yjpj  +  cr2,  Im  r  =  err, 
j= i 

so  that  Re  r  is  even  and  Im  r  is  odd  on  the  image  contour  T'.  Note  also  that 

dr  _  ^  hjp 

dp  J=i  y/pij  -  p  2 


(0  ?0) 


In  particular, 


dr 

dp 


=  r  >  0. 


(0.27) 


(0.28) 


p=0 


From  this  it  follows  that  arg  dr  —  7r/2,  because  arg  dp  =  7r/2  at  p  =  0.  Thus,  the 
contour  is  directed  vertically  upward  at  the  image  point 


t(0)  =  Y,  hiPh 

j= i 


(0.29) 


which  is  the  travel  time  along  a  vertical  raypath  from  the  source  depth  to  the  obser¬ 
vation  depth.  That  path  includes  all  multiple  reflections  at  interfaces  that  are  implied 
by  the  terms  of  the  sum  in  w(p,r,z). 

In  general,  on  this  image  path 


Re  r  — »  |<j  I  Y,hJ’  Im  r  — >  or,  \a\  — ♦  oo, 
j=  i 


and 


Im  r 


sgn(fr) 


\a\  — *  oo. 


(0.30) 


(0.31) 


Rer  °  '  '  £/=i  hj  ’ 

That  is,  the  image  contour  has  an  asymptote  in  the  r-plane  with  slope  proportional 

to  the  ratio  of  the  horizontal  travel  distance  to  the  vertical  offset.  In  the  simplest 
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case,  J  =  1,  this  curve  is  one  branch  of  a  hyperbola.  For  larger  values  of  J ,  the  basic 
hyperbola-like  nature  of  a  symmetric  curve  with  symmetric  asymptotes  still  holds. 
Figure  0.2  is  an  example  of  an  image  contour  T'  for  a  case  in  which  J  =  2  and  the  two 
propagation  speeds  are  3km/sec  and  4km/s  with  h\  =  h?  =  1.2km  and  r  =  2.5km. 

Because  the  classic  Cagniard  integration  path  in  the  complex  p-plane  is  also 
hyperbolic-like,  the  reader  is  cautioned  not  to  confuse  it  with  the  previously  dis¬ 
cussed  image  contour  T'  of  the  imaginary  p- axis.  This  path  is  in  the  r-domain;  the 
Cagniard  path  is  in  the  p-domain. 


The  Saddle  Points  of  w 


Saddle  points  are  defined  as  those  points  at  which  dw/dp  =  dr /dp  =  0.  From 
(0.27),  one  can  verify  that  there  is  always  a  point  on  the  Rep  axis  where  this  is  true. 
In  that  equation,  one  can  see  that 


— - »  — oo,  as  p  — ►  min(pj). 

dp  j 


(0.32) 


Because  dr/dp  is  positive  at  p  =  0,  from  (0.28),  dr/dp  must  pass  through  zero 
somewhere  on  this  interval.  This  is  the  saddle  points  to  be  denoted  by  p3.  At  this 
point, 


r 


j 


=  E/ti 

i= i 


(0.33) 


This  equation  defines  the  geometrical  optics  ray  that  connects  the  origin  with  the 
point  (x,  z)  including  multiple  reflections  and  refractions  satisfying  Snell’s  law.  To 
see  this  latter  point,  observe  that  Snell’s  law  for  reflection  and  refraction  guarantees 
that 


pa=Pjsindj,  j  =  1,  J, 


(0.34) 


is  the  same  for  each  j  (even  for  mode  converted  waves).  Equation  (0.32)  assures  that 
p.,  <  min j(pj):  and  thus,  that  each  6j  is  real.  The  horizontal  travel  in  the  jth  layer  is 
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given  by  hj  tan  Oy,  thus,  the  total  horizontal  travel  is 

j 

r  =  hj  tan  Qj.  (0.35) 

i= i 

If  we  use  (0.34)  to  define  tan0_,  in  terms  of  ps  and  pj ,  the  result  is  (0.33). 


The  saddle  point  ps  is  a  significant  point  for  our  analysis.  At  this  point,  the 
mapping  from  p  to  r  ceases  to  be  conformal;  that  is,  the  angle  between  two  intersecting 
arcs  in  the  p-domain  will  not  be  preserved  in  the  r-domain.  To  determine  how  this 
angle  is  changed,  it  is  necessary  to  find  the  first  nonvanishing  derivative  at  ps.  To 
this  end,  we  differentiate  (0.27)  to  obtain 


d2', 


=  -E 


hjpj 


dp2  fr[  [p2  _  p2j 


3/2’ 


(0.36) 


The  explicit  formula  shows  that  d2T/dp2  is  negative  on  the  interval  from  0  to  miiij(p_,). 
This  establishes  that  p  —  ps  is  the  only  saddle  point  on  this  interval  and  that  the 
second  derivative  is  negative  at  this  saddle  point.  Let  us  set 


r{ps)  ~  ts. 


(0.37) 


Because  d2r/dp 2  is  nonzero, 


r  -  r. 


1  d2r(p3 ) 

2  dp2 


{P-P*)2< 


(0  38) 


from  which  it  follows  that  changes  in  argument  of  p  —  p9  are  doubled  in  the  r-domain. 


Of  particular  interest  to  us  are  the  paths  through  ps  on  which  Imre  =  constant. 
In  on-  case,  that  constant  is  zero,  since  w(ps)  is  real.  Of  necessity,  the  image  of  these 
paths  must  be  a  segment  of  the  Rer  axis,  and  these  paths  Imre  =  0  through  the 
saddle  point  are  segments  of  the  Cagniard  path. 

In  the  method  of  steepest  descents,  these  paths  of  constant  Im  w  are  called  paths 
of  steepest  ascent  or  descent.  The  reason  is  that  when  we  think  in  terms  of  the  surface 
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Re  w,  these  are  the  paths  of  maximal  change  in  Re  w  away  from  the  reference  point,  or 
literally  the  directions  of  steepest  ascent  and  descent  on  this  surface.  In  the  method 
of  steepest  descents,  the  paths  of  steepest  descent  are  the  paths  of  interest  (since  the 
exponent  of  interest  is  +sw):  here,  the  paths  of  choice  are  the  paths  of  steepest  ascent 
(since  the  exponent  of  interest  is  —  sw)  because  in  either  case  we  are  interested  in  the 
paths  on  which  the  integrand  decays  exponentially  to  zero  at  a  maximal  rate.  We 
will  denote  by  S  the  contour  made  up  of  the  two  paths  of  steepest  ascent  through  ps 
oriented  in  the  direction  of  ascent  in  the  upper  half-plane.  Figure  0.2  shows  S'  for 
the  same  example  considered  in  Figure  0.1.  Note  that  Figure  0.2  shows  a  contour 
in  the  r-domain,  while  Figure  0.1  shows  a  contour  in  the  p-domain.  The  contour 
T'  in  Figure  0.2  is  the  image  of  the  Imp  axis.  As  noted  above,  the  image  S'  of  the 
contour  S  in  Figure  0.1  is  the  doubly  covered  Rer  axis,  starting  at  r,  and  extending 
to  infinity.  This  will  be  verified  below.  When  viewed  as  a  conformal  mapping,  the 
image  of  the  region  between  the  T  and  S  in  the  p-domain  (Figure  0.1)  is  the  region 
between  T'  and  S'  in  the  r-domain  (Figure  0.2).  As  noted  in  the  introduction,  we 
refer  to  this  former  region  of  the  p-plane  as  the  primary  p-domain  and  we  refer  to  its 
image  in  the  r-domain  as  the  primary  r-domain. 

On  5,  Imre  =  0,  as  required;  however,  there  are  two  curves  through  ps  with  this 
property  and  they  are  orthogonal  to  one  another.  How  can  we  see  that  the  image  of 
the  contour  indicated  here  is  S',  as  claimed?  Note  first,  from  (0.36)  and  (0.38)  that 

rs  -  r  as  A(pa  -  pf  (0.39) 

with  A  >  0.  Here,  reversing  the  order  of  r's  on  the  left  side  of  the  equation  has 
allowed  us  to  express  the  approximation  (0.38)  in  terms  of  a  positive  constant,  A. 
Reversing  the  order  of  the  p's  on  the  right  side  of  the  equation  allows  us  to  express 
the  right  side  in  terms  of  an  expression  for  which  we  know  the  argument  at  p  =  0  and 
then,  by  continuity  allows  us  to  know  the  argument  everywhere.  In  particular,  one 
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can  verify  that  when  p  moves  from  the  origin  to  a  point  on  S  above  ps  and  nearby, 
arg  (ps  —  p)  «  -  -  7T  /  2  and  arg(rs  —t)zz  —tt.  Similarly,  when  p  moves  from  the  origin  to 
a  point  on  S  below  ps  and  nearby,  arg(p„  —  p)  «  +7r/2  and  arg(r,  -  r)  «  +7T.  In  either 
case,  ts  —  t  is  negative  as  required.  One  can  further  verify  that  on  the  alternative 
paths  of  constant  Im  u>,  Re  w  decreases.  That  is,  the  image  of  these  paths  is  the  Rer 
axis  to  the  left  of  ts. 


Of  special  interest  is  the  Rep  axis  to  the  right  of  ps.  Here  we  must  have  arg(p5  — 
p)  =  ±7r,  depending  upon  whether  this  region  was  reached  by  p  passing  above  (  — ) 
or  below  (  +  )  ps  when  starting  from  p  =  0.  It  follows  from  (0.39)  that  for  these  cases 
arg(r«  —  r)  =  ±27r,  which  are  choices  of  the  argument  outside  of  the  primary  r- 
domain.  These  values  of  r  correspond  to  points  on  the  second  Riemann  sheet  of  the 
mapping  from  p  to  r  and  values,  t  <  ts. 


It  will  prove  useful  to  evaluate  ts.  To  this  end,  we  use  (0.33)  in  (0.23)  to  find 


Ts  =  E  k3 


P2s 


^  3  n — 2 

yjPj  -  P2s 


+  hjy/p'j  ~  P2 


(0.40) 


=  E  }h 


p) 


U  VpF^' 


Let  us  define 


rj  =  -fhfS  Pi  =  \lr]  +  hy  (0-41) 

y/Pj  ~Ps 

That  is,  r j  is  the  segment  of  the  offset  r  at  p  =  ps  in  the  j-tli  layer  and  pj  is  the 
length  of  the  ray  segment  in  the  j-th  layer.  Substitution  of  the  definition  of  rj  into 
the  definition  of  Pj  yields 

/jjT).- 

(0.42) 


_  hjPj 
Pj  — 


\/p2-P2' 


and 


rs  —  XI  PjPj 

J=I 


(0.43) 
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is  i,he  total  traveltime  along  the  ray.  When  we  compare  this  result  with  (0.29),  we 
see  that 

r,  >  r( 0),  (0.44) 

with  equality  holding  only  for  r  =  0.  Indeed,  from  (0.36),  r,  is  the  maximum  of  r  for 
0  <  p  <  min j(pj)  and  this  inequality  holds  for  r(0),  replaced  by  r(p)  for  any  p  in  this 
interval. 


The  point  ts  is  a  branch  point  of  the  integrand  in  (0.3).  To  see  why  this  is  so,  let 
us  solve  (  0.38)  for  p: 

2(r  —  Ts) 

P-P*  ~  \  -3TTTTT-  (°-4°) 


\J  d*rs/dp 2  ‘ 


Differentiation  of  this  equation  yields 

dp  1 

*  cv  _ _ 

dT  yj2{r  -  Ts)(pTs/dp 2 


(0.46) 


which  verifies  the  claim. 


We  can  draw  an  important  conclusion  from  this  observation.  The  image  in  the 
T-domain  of  the  steepest  descent  path  in  the  p-domain  is  a  segment  of  the  Rer  axis 
that  lies  to  the  right  of  the  image  of  the  imaginary  axis  as  shown  in  Figure  0.2. 

In  the  T-domain,  we  are  most  concerned  with  the  region  in  which  0  <  arg(r  — r5)  < 
2n  that  lies  in  the  primary  T-domain.  To  see  this,  note  from  (0.36)  that  dPr/dp2  <  0 
and,  of  necessity,  argp  —  ps  =  7r/2,  37t/2,  for  t  —  ts  to  be  positive.  When  we  take 
arg  cPr/dp2  =  —  7r,  we  find  that  these  two  directions  correspond  to  argT  —  rs  =  0,  2n, 
respectively.  Then,  in  this  neighborhood,  values  of  a.Tgp  —  ps  less  than  7r/2  or  greater 
than  37t/2  must  correspond  to  images  points  with  argT  —  t,  less  than  zero  or  greater 
than  27r,  respectively.  These  regions  correspond  to  another  Riemann  sheet  of  the 
multi-sheeted  surface  of  the  T-domain. 


Of  particular  importance,  singular  points  of  f(p)  in  the  primary  p-domain  map 
into  points  in  the  primary  T-domain.  As  previously  noted,  not  all  of  these  need  to 
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be  singular  points  of  the  integrand  f{r)dpjdr  in  (0.3).  For  those  that  are  singular 
points  after  transformation,  deformation  of  the  contour  T'  onto  the  Rer  axis  cannot 
proceed  without  including  loop  contours  around  these  points.  On  the  other  hand, 
singularities  to  the  right  of  the  path  of  steepest  descent  in  the  p-domain  and  out  of 
the  primary  p-domain,  must  of  necessity  map  onto  these  other  sheets  of  the  Riemann 
surface. 

The  Branch  Points  of  w 

We  turn  now  to  consideration  of  the  points  pj  which  are  the  branch  points  of  w(p) 
in  (0.2).  In  the  neighborhood  of  pj, 

T  -  r(pj)  «  hjyjlpjyjpj  -  p  +  0(j)j  -  p),  (0.47) 

from  which  it  follows  that 

P i  -  P  «  (r-^(ff>)2  +  0((T  -  r(Pi))3).  (0.48) 

From  the  first  ecpiation  here,  we  conclude  that  y/pj  —  p  is  linear  in  r  —  r(pj).  Thus, 
if  this  square  root  appears  in  /(p),  it  is  not  a  branch  point  of  /(p(r)).  From  the 
second  equation  here,  we  conclude  that  dp/dr  is  also  linear  in  r  —  r(pj).  Thus,  the 
branch  points  of  iv(p)  are  not  a  priori  branch  points  of  the  transformed  integrand 
f{p(r))dp/dT  of  (0.3). 

On  the  other  hand,  other  singular  points  of  /(p)  must  of  necessity  be  singularities 
of  corresponding  type — branch  point,  pole,  etc. — of  f(p(r)),  because  drjdp  is  finite 
and  nonzero  at  such  points. 

The  Choice  o  =  tt/2. 

We  close  this  section  by  returning  to  the  question  of  the  choice  of  final  argument 
tt/2  for  w.  For  any  other  choice,  the  path  F  would  not  be  vertical.  Then,  for  r  near 
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enough  to  zero,  F  would  intersect  5,  in  which  case  the  image  domain  in  r  becomes 
much  more  complicated.  In  particular,  the  image  path  T'  would  cut  through  the 
image  path  S'  and  part  of  the  primary  r-domain  would  be  on  the  second  Riemann 
sheet. 

On  the  other  hand,  it  has  been  shown  that  knowing  ihe  Laplace  transform  for 
real  values  of  s  suffices  for  the  determination  of  its  inverse  function  (technically,  up 
to  a  function  of  measure  zero,  but  uniquely  at  points  where  the  inverse  transform 
is  at  least  left-  or  right-continuous).  This  follows  from  the  observation  that  the 
Laplace  transform  must  be  analytic  in  some  right-half  s  —  plane  and  Lerch’s  theorem 
[Sneddon,  1972,  Widder,  1959],  which  assures  uniqueness  of  the  Laplace  inversion  of 
such  functions.  Then,  one  need  only  observe  that  an  analytic  function  is  determined 
in  its  domain  of  analyticity  by  its  values  on  a  line  segment  in  that  domain,  no  matter 
how  small.  Any  segment  of  the  Res  axis  will  do.  It  is  standard  in  the  literature 
to  simply  state  that  knowledge  of  the  Laplace  transform  on  the  real  axis  suffices  to 
determine  the  inverse  function  (in  the  time  domain)  uniquely. 

EXAMPLES 

Here  we  will  demonstrate  the  application  of  the  Cagniard-de  Hoop  method  on  var¬ 
ious  examples  designed  to  demonstrate  the  analysis  above  and  bring  out  the  features 
of  this  interplay  of  two  complex  planes.  The  reader  should  note  that  the  image  of 
the  steepest  ascent  path  in  the  /^-domain  is  always  at  least  a  part  of  the  loop  integral 
path  around  the  real  axis  in  the  r-domain  and,  sometimes,  may  be  all  of  it.  Whether 
or  not  it  is  all  of  that  path  depends  on  whether  or  not  there  are  other  singularities  of 
/(/>)  in  the  primary  p-domain  other  than  branch  points  of  order  one-half  at  the  same 
locations  at  the  branch  points  of  the  exponent  w(p,r,z). 
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Two-dimensional  free  space  Green’s  function 


In  this  first  example,  the  only  singularity  of  the  integrand  in  the  r-domain  will  be 
the  branch  point  at  r,.  We  will  verify  further  that  the  integrand  is  analytic  at  the 
image  of  the  (single)  branch  point  of  w(r,z,p). 


We  consider 


(0.49) 


This  is  an  integral  representation  of  the  two-dimensional  Green's  function  for  the 
Helmholtz  equation  for  acoustic  waves  in  a  constant  density  medium. 


V'-\,  +  =  -b{x)b{z). 

(T 


(0.50) 


Here  V*  is  the  two  dimensional  Laplacian.  See  Bleistein  (1984),  for  example. 


This  representation  is  a  special  case  of  (0.6)  and  the  transformation  of  the  solution 
proceeds  as  in  the  discussion  below  (6)  to  yield  (0.1)  with  w  defined  by  (0.2)  and 
.7  =  1.  That  is. 

s)  =  — exp  {-snip.  |.r|,r)}  (0.51) 

4“'  J]  \Jpi~P 2 

with  ir  given  by 

«'(/>•  |-r|. -)  =  }>\-r\  +  yjp\  —  p2|r|,  p,  =  1/(4.  (0.52) 


In  this  examj)le.  the  only  critical  points  of  the  integrand  are  the  branch  points  at 
±p\  and  the  saddle  point  determined  by  (0.33)  with  r  =  |.r|: 

p,  =  y/,1  sinflj.  .r  =  psin#,  c  =  pcosW,  p  =  v t*  +  .  (0.53) 

We  need  no  index  on  ft  here  because  there  is  only  one  value,  pt . 

The  transformation  (0.23)  in  this  case  is 

r  =  »-(p.  |.r|,;)  =  p|.r|  4-  yjp\  —  p2|~|  (0.54) 
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and  (0.24)  becomes 


1  [  dp  dr 

[r,  c,  s)  -  — :  /  - - 7= - exp{-ST}, 

4rr/  Jr'  dr  pi  _  ^2 


(0.55) 


with  T'  as  shown  in  Figure  0.2  except  that  the  specific  numbers  of  that  example  are 
not  relevant  here.  The  theory  assures  us  that  the  integrand  has  no  singularities  except 
for  the  branch  point  at 

ts  =  pip,  (0.56) 

and  that  there  arc  no  other  singularities  of  the  integrand  in  (0.24).  We  will  explicitly 
verify  that  below. 


We  can  now  close  the  contour  of  integration  T'  around  the  branch  point  at  r, 
to  obtain  a  loop  integral  on  the  contour  S'  as  shown  in  Figure  0.1.  We  denote  the 
solution  p(r)  for  r  on  the  upper  side  of  the  loop  by  p+  and  the  solution  p(r)  for  r  on 
the  lower  side  of  tin1  loop  by  p_.  Then  we  can  write 


u(r.z.s)  = 


47 ri 


dp+  1 

_,lT  vVf - p2+ 


fip_  1 

dT  \[p{  ~  P-. 


exp  {  — sr}t/r. 


(0.57) 


Now  the  time-domain  solution  can  be  read  off  as  the  integrand  here.  However,  there 
is  one  further  simplification  that  comes  from  the  fact  that  the  two  terms  in  the 
integrand  are  complex  conjugates  of  each  other.  This  fact  observation  follows  from 
the  Schwarz  reflection  principle  [Levinson  and  Redheffer,  1970],  which  assures  us  that 
an  analytic  funct  ion  must  assume  complex  conjugate  values  across  a  line  in  its  domain 
of  analyticifv  on  which  the  function  is  real.  Here  the  line  along  which  the  function 
p(  r)  is  real,  is  the  line  ;n|~|  <  r  <  r,  =  pip,  which  is  the  image  of  0  <  p  <  p„  where 
we  know  that  ir  is  real.  Thus,  complex  conjugate  points  in  either  domain  correspond 
to  complex  conjugate  points  in  the  other  and  we  conclude  that 


l'(jr.z.t)  =  Im 


<h>+  1 

<lT  2?r yjp\  -  pi 


H(t-r ,)■ 


(0.58) 
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In  this  equation,  H(t)  is  the  Heaviside  function ,  equal  to  unity  for  positive  argument 
and  zero  for  negative  argument. 

For  this  particular  case,  we  can  explicitly  solve  fo rp(r)  because  (0.54)  is  equivalent 
to  a  quadratic  equation.  The  solution  is 


P  = 


k|r  -  \~\\Ip\p2  ~  t2  Ip2- 


(0.59) 


In  this  solution,  it  was  necessary  to  make  a  choice  of  sign  of  the  square  root.  We 
chose  the  sign  for  which 

P  =  0  <=>  r  =pi\z\, 


which  is  consistent  with  (0.54).  The  other  choice  of  sign  of  the  square  root  corresponds 
to  the  mapping  of  p  to  the  second  Riemann  sheet  of  the  r-domain.  We  see  here  that 
the  only  singularities  of  this  function  are  at  r  =  ±pip  =  ±ts.  One  of  these  is  as 
predicted  by  the  theory  and  the  other  is  completely  outside  the  domain  of  interest, 
since  it  is  not  in  the  primary  r-domain. 

One  can  further  confirm  that 


\z\t  +  \x\jp\p2  -  T2 


(0.60) 


which  also  is  singular  only  at  ±rs  and  not  at  r(pl)  =  pi|r|.  This  is  as  predicted  by 
the  theory. 

Finally, 

^  =  +  (061) 

P2\Jp\p2  ~  T2 

Since  the  Heaviside  function  in  (0.58)  confines  the  nonzero  portion  of  the  solution 
to  t  >  t,  =  ]>ip,  and  since  we  wrote  this  result  in  terms  of  p+,  we  must  determine 
what  happens  to  the  square  root  in  (0.59)  when  we  pass  over  the  branch  point  at 
p\p  in  the  r-domain.  What  is  crucial  here  is  that  arg(pi/r  —  r)  varies  from  0  to  — 7r 
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on  this  trajectory.  We  conclude  then  that  the  square  root  must  have  an  argument  of 
—  7t/2  (essentially,  a  multiplier  of  —i  on  a  real  positive  square  root).  That  is, 


P+  = 


x\t  +  i\z\\J  T2  —  p\p2 


Ip2,  t  >  pip. 


(0.62) 


However,  we  need  not  deal  separately  with  the  derivative  and  /(p+(r))  because  it 
follows  from  (0.60)  and  (0.61)  that 

dp+  1  i 


dr 


y/rf-  P+  y/r2  -p\p2 


,  T  >  Pip. 


(0.63) 


iow  we  can  evaluate  (0.58)  explicitly  to  obtain 

U(x,z,t)  = 


(0.64) 


2n^t2  —  p2/c\ 

Here  we  return  to  the  use  of  the  propagation  speed  c\,  rather  than  the  slowness  p\ . 
This  is  the  two-dimensional  free  space  Green’s  function  for  the  wave  equation  (0.50). 

Point  Source  over  a  Half  Space 

As  a  second  example  we  consider 

1  r°°  dk  / - 

v(x,z,u)  =  —  - — :  /  —j===:R(k/u;)exp{i[kx  +  iJu)2/c2  —  k2(2h  —  z)]}, 

4tti  J-oo  yju)2/c2i  -  k2 

(0.65) 


with 


R(k/u) 


_  c /cj  -  viu?  -  yi/ci  - 


(0.66) 


k2/u>2  +  ^jl/c2  —  k2/u 2 
This  is  the  upward  scattered  acoustic  wave  in  a  constant  density  medium  for  a  point 

source  at  height  h,  over  a  half  space.  See  Bleistein  (1984),  for  example.  This  solution, 
when  added  to  the  previous  one,  provides  the  total  response  above  the  interface 
between  two  homogeneous  half  spaces  with  different  acoustic  speeds. 

This  solution  again  fits  the  form  of  (0.6)  and  the  transformation  of  the  solution 
proceeds  jus  in  the  discussion  below  that  equation  to  yield  (0.1)  with  w  defined  by 
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(0.2)  and  J  —  1.  That  is, 


with 


u(x,z,s)  -  f  -jJ==R(p)exp{-sw(p,\x\,z)}, 
47H  Jr  Y'PX  —  p2 


w(p,  Ml  *)  =  pM  +  ]/pf  -p2(2h  -  z), 

jhp) = 

yj?\  ~P2  +  \jp\  ~  P2 
Pi  =  1/ci,  p2  =  l/c2. 


(0.67) 


(0.68) 


The  new  feature  of  this  example  is  that  /(p)  now  has  a  branch  point  at  P2,  which 
is  not  a  branch  point  of  the  exponent  w.  As  in  the  previous  example,  the  branch 
point  at  pi  appears  to  the  same  order  in  /(p)  and  w(p ,  |.t|,  z). 

Note  that,  except  for  replacing  j.r|  by  2 h  —  z,  the  analysis  of  the  transformation 
from  p  to  t  is  exactly  as  in  the  previous  example.  In  particular,  the  saddle  point 
determined  by  (0.33)  with  r  =  |x|  becomes 

Ps  =  Pi |  sin0|,  x  =  p sin0,  2h  —  z  =  pcosd,  p  =  yjx2  +  (2h  —  z)2.  (0.69) 

The  region  in  which  2 h  —  z  >  0  corresponds  to  \6\  <  7r/2. 

The  transformed  integral  (0.24)  for  this  case  is 

"(-r’'’'s)  =  hL  cxp(‘ST)’  (0-70) 

with  T'  as  shown  in  Figure  0.2,  except  that  again,  the  specific  numbers  of  the  figure 
are  not  relevant  here. 

Wo  must  be  concerned  with  the  location  of  the  branch  point  p2  and  its  image 
in  the  r-domain.  In  particular,  we  must  know  when  the  image  falls  in  the  primary 
r-domain.  If  p2  >  p\,  this  cannot  happen,  since  p.,  <  p\  =  min{pi,p2);  see  preceding 
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Fig.  0.3.  Physical  domain  showing  pre-  and  post-critical  regions  defined  in  terms  of 

e. 

discussion  on  the  location  of  the  saddle  point.  Therefore,  the  interesting  case  occurs 
when 

P2  <  Pi  or  ci  <  c2.  (0-71) 

Of  course,  the  reader  familiar  with  this  problem  will  observe  that  this  is,  indeed,  the 
case  in  which  head  waves  arise  and  that  they  manifest  themselves  in  the  solution 
through  the  effect  of  this  second  branch  point.  We  proceed  with  this  more  interesting 
case. 

We  introduce  the  critical  angle,  0C  by  the  equation 

P2  =  Pisin0c  (0.72) 

and  note  that  p2  will  be  in  the  primary  p-domain  if 

sin0c  <  |  sin0|.  (0.73) 

From  Figure  0.3,  we  can  see  that  this  is  the  post-critical  reflection  region  in  which 
the  incident  wave  is  totally  reflected.  Therefore,  for  this  physical  region  we  now  know 
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Fig.  0.4.  Deformed  path  of  integration  in  the  r-domain  for  p 2  in  the  primary  p- 
domain. 


that  the  image  point  rc,  given  by 


tc  =  iu(p2)  =  ppi{sin#c|  sin#|  +  cos0ccos0} 


=  PPicos{0c  -\e\)  =  (p/cO  cos{0c  -  |0|}, 


(0.74) 


lies  in  the  primary  r-domain.  Thus,  the  deformation  of  T'  in  this  case,  must  be  a 
loop  integral  around  rc  as  shown  in  Figure  0.4. 

To  see  the  advantage  of  using  two  complex  planes  in  the  Cagniard-de  Hoop 
method,  we  note,  as  in  the  preceding  example,  that 


T  =  r{p.,)  =  w(p3,  |a:|,  z)  =  p/Cl  >  rc  =  (p/ci){cos0c  -  |^|}  =  rc,  (0.75) 

and  this  identity  holds  no  matter  whether  P2  is  in  the  primary  p-domain  or  not. 
Again,  note  that  r,  is  the  maximum  of  r  for  0  <  p  <  p\.  Therefore,  by  just  knowing 
the  value  of  rc  we  cannot  tell  whether  or  not  it  lies  between  T'  and  S'.  It  is  only  when 
we  sec  that  it  is  the  image  of  a  point  in  the  primary  p-domain  that  we  can  be  sure 
that  it  lies  between  these  image  contours.  When  p2  is  to  the  right  of  5,  wc  know  that 
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its  image  in  the  r-domain  lies  on  the  second  Riemann  sheet  of  the  branch  point  at 
ts  and  not  between  the  image  contours,  H  and  S'.  However,  by  understanding  the 
mapping  from  p  to  r  we  readily  know  that  for  P2  to  the  right  of  the  primary  p-domain, 
its  image  must  be  on  the  lower  Riemann  sheet  with  respect  to  the  branch  point  at  rs 
and  at  some  r  <  ts.  Thus,  such  a  choice  of  P2  does  not  interfere  with  the  deformation 
of  R  onto  S1 . 

We  now  have  the  machinery  in  place  to  find  a  time  domain  representation  of  the 
upward  propagating  wave  for  this  problem  in  both  the  pre-  and  post-critical  regions 
in  Figure  0.3.  We  \. ill  discuss  them  in  that  order. 


Case  I:  \9\  <  6C. 


From  (0.69)  and  (0.72),  it  follows  that  in  this  case  ps  <  p2  and  the  latter  point  is 
outside  the  primary  p-domain.  Therefore,  its  image  is  outside  the  primary  r-domain 
and  the  deformation  from  T'  to  S'  proceeds  unimpeded  by  the  presence  of  the  image 
rc  of  p2-  In  analogy  with  the  transformation  from  (0.51)  to  (0.57),  in  this  case  (0.67) 
becomes 


u(x,  z,  s ) 


1  r°° 

dp+ 

R(p+) 

dp _ 

Rip-) 

47 xi  Jt, 

dr  , 

Jv\  -  P+ 

dr 

\Jp\~p1-. 

exp  {—srjdr 


(0.76) 


and 

U(x,z,t)  =  Ira  7°=  ,  \mt-r,).  (0.77) 

l  dT  2nyjpi  -  p2+  J 

For  this  case,  (0.63)  is  still  valid,  except  that  p  is  now  defined  by  (0.69).  Since  the 
right  side  of  this  last  expression  is  purely  imaginary,  we  extract  the  imaginary  part 
in  (0.77)  by  choosing  the  real  part  of  R(p+ );  that  is, 


U(x,z,t)=  H{)  />/C'bRe  (B(p+)).  (0.78) 

2-KsJt2  -  p2lc{ 

This  result  agrees  with  Aki  and  Richards,  I,  1980,  p.  230,  eq.  (6.57). 
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Case  II:  6  >  6C. 


In  this  case,  as  previously  noted,  p2  <  ps  and  lies  in  the  primary  p-domain,  its 
image  rc  lies  in  the  primary  r-domain  and  the  deformed  contour  in  the  7-domain  is 
as  shown  in  Figure  0.4.  Consequently,  we  replace  (0.76)  by 


u(x,z,s) 


1 

47ri 


r + r 

dp+ 

R{P+) 

dp_ 

R(p- )  ' 

Jtc  Jtb 

dr 

\/pi  ~  P+ 

dr 

Jpi-pl 

exp  {—sr}dT 


(0.79) 


=  «i(x,2,s)  +  u2(x,z,s), 


with  u i  and  u2  corresponding  to  the  first  and  second  integrals,  respectively. 

Let  us  first  consider  the  integral  ui(x,z,s)  and  its  inverse  transform  Ui(x,z,t). 
The  fact  that  both  limits  of  integration  are  finite  tells  us  that  the  inverse  transform 
is  a  function  that  “turns  on”  at  t  =  tc  and  “turns  off”  at  t  —  rs.  This  finite  interval 
of  integration  is  the  image  of  the  line  segment,  p2  <  p  <  ps  in  the  primary  p-domain. 
On  this  interval,  the  change  of  variables  defined  in  (0.68)  is  real.  In  analogy  with 
(0.60)  and  (0.61),  we  find  here  that 


\Jpi  -p2  =  [2/t  -  z\t  -I-  \x\\Jp\p2  -  t 2 


Ip 2, 


(0.80) 


and 


dp  _  [2/t  -  z]t  +  \x\sjp\p2  -  t2 
dr 


P2\]p2\P 2  -  t2 


with  p  defined  by  (0.69).  Thus,  we  conclude  that 


(0.81) 


=  Im  <  %— ?===  \[H(1  -  t,)  -  H(t  -  r,)] 


dr 


2^ 


pi-  p2 


Im  {R(p+)}  Hit  —  tc)  —  H(t  —  ts) 

277  sjp2/ci  -  t 2 


(0.82) 
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For  u2(x,  z,  s ),  the  path  of  integration  is  one  for  which  the  p(r)  is  again  complex 
valued  with  the  derivative  as  in  the  above  case  9  <  6C, 

U2(x,z,t)  =  -/M  Re{i?(Pt)}.  (0.83) 

2't/l  -  P2/c2, 

These  results  agree  with  Aki  and  Richards,  I,  p.  234,  eq.  (6.58). 

This  completes  the  discussion  of  this  elementary  problem.  It  is  typical  of  the  type 
of  analysis  that  is  necessary  when  viewing  the  change  of  variables  as  a  mapping  from 
one  complex  plane  to  another. 


Three-dimensional  free  space  Green’s  function 


We  now  consider 


1  roc  roc  dkidk2exp{i[kix  +  k2y  +  Ju)2/c^  —  k*  —  kl\z\)} 
v(x,y,z,u)  = /  - - - - 

oil “l  J—oo  7—00 


yju)2/^  -  -  k2 


(0.84) 


This  is  the  three-dimensional  analog  of  (0.49)  and  is  of  the  form,  (0.14),  with  J  =  1 
and 

11  11 


g{k  1,  k2,u)  —  — ; 


87r2*  \J oS2  j c\  -  k\  -  kl  8n2luJ  \J  1  /cj  -  k\/uj2  -  k2 / u2 
so  that  in  (0.15),  n  =  —  1  . 


,  (0.85) 


The  transformations  below  (0.15)  can  be  applied  here,  leading  to  functions  /  and 
v  defined  by  (0.21).  That  is, 

11  11 


lip)  = 


87r2*'  ^/l/c2  +  q2-p2  \]p\  ~  P2 


(0.86) 


and 


v(x,  y,  2,  is)  =  su(r,  2,  s).  (0.87) 

In  the  first  equation  pi  is  defined  by  (0.22)  and  in  the  second  equation,  multiplication 
by  s  means  that  we  must  differentiate  with  respect  to  t  in  addition  to  integrating 
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with  respect  to  q  in  order  to  obtain  the  time  domain  Green’s  function  from  the  result 
of  applying  the  Cagniard  method  to  the  function  u. 

With  these  results  we  find  that  we  must  analyze  the  integral  representation 


"(w)  =  iklq^exp{-sw{p'r'2)h 


(0.88) 


with  w  given  by 


w(p,  r,  z )  =  pr  +  yjp\  -  p2\z\,  p\  =  1  / c\  +  q 2. 


(0.89) 


These  equations  are  the  same  as  (0.51)  and  (0.52)  except  for  (i)  replacement  of  x  and 
|x|  by  r  and  (ii)  an  extra  factor  of  1/27T  in  the  definition  of  u.  Thus,  we  read  off  the 


result  of  the  Cagniard  method  from  (0.64)  as 


w(t-pAi) 


IT!  r  tl  2 

U(r,z,t)  =  —  ^  p  =  vr2  +  *. 


(0.90) 


Here  we  see  the  power  of  de  Hoop’s  transformation:  except  for  the  quadrature  and 
differentiation  to  follow,  the  application  of  the  Cagniard  method  to  the  represen¬ 
tation  of  the  three  dimensional  Green’s  function  is  reduced  to  the  analysis  of  the 
two-dimensional  problem  treated  earlier. 

As  noted  above,  it  is  necessary  to  integrate  this  result  with  respect  to  q  : 


LU{r'-'t)dq  =  4A,  -  J  _ 


I It2  -  P2q2 


(0.91) 


H{t-  p/ci) 
47T  p 


Finally,  we  differentiate  with  respect  to  t  to  obtain  v  : 


v(x,y,z,t)  = 


6(t  -  p/cx) 


(0.92) 
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In  Chapter  5  of  Cagniard’s  (1939,  1962)  book  and  in  Dix  (1954),  polar  coordinates 
are  used  in  the  representation  of  the  Green’s  function.  This  leads  to  a  pair  of  integrals 
in  transverse  wave  number  and  angle  between  the  transverse  wave  vector  and  ( x,y ). 
The  first  of  these  integrals  is  treated  bv  Cagniard’s  method.  However,  the  function  w 
depends  on  the  polar  angle,  making  the  subsequent  integration  (the  analog  of  the  q 
integral  here)  much  more  difficult.  This  latter  integration  is  a  tour  de  force  in  complex 
function  theory  by  those  authors.  The  de  Hoop  transformations  avoid  this  extremely 
difficult  analysis,  making  application  of  his  modification  of  Cagniard’s  method  much 
more  accessible  and  a  method  of  choice  for  a  broad  class  of  problems. 

The  Sommerfeld  Half-Plane  Problem 

We  will  consider  here  diffraction  of  a  plane  wave  by  a  half-plane  in  the  case  where 
there  is  no  equivalent  causa!  problem.  This  is  not  a  problem  which  is  normally  treated 
by  the  Cagniard-de  Hoop  method  starting  from  the  Laplace  transform,  because  that 
method  requires  an  equivalent  causal  problem  for  the  (one-sided)  Laplace  transform 
to  make  sense.  However,  the  structure  of  the  problem  certainly  lends  itself  to  treat¬ 
ment  by  the  Cagniard-de  Hoop  method,  if  we  start  from  a  Fourier  representation  of 
the  solution  and  carefully  observe  the  interplay  of  singularities  with  the  rotation  of 
contour  associated  with  the  transformation  from  {k,u)  to  (p,  s). 

The  underlying  problem  is  as  follows.  A  plane  acoustic  wave  in  a  constant  density 
medium  is  incident  from  infinity  on  a  half-plane  occupying  the  region,  x  >  0,  2  =  0. 
The  wave  is  also  parallel  to  the  y  axis  so  that  the  problem  is  two-dimensional.  Given 
a  boundary  condition,  the  objective  is  to  find  the  total  field  everywhere.  We  will 
assume  the  boundary  condition  that  the  total  field  is  zero  on  the  scatterer. 

In  Figure  0.5,  we  show  two  possible  directions  of  incident  wave  from  above.  The 
waves  are  distinguished  by  the  relative  inclination  of  the  incidence  direction,  as  defined 
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Fig.  0.5.  Plane  wave  incident  on  a  half-plane;  two  possible  cases  shown,  distinguished 
by  relative  inclination  to  the  scatterer.  The  angle  of  incidence  is  n  —  9C. 

by  the  angle  z  -  6C.  For  0C  >  tt/2 — incidence  from  the  left — one  can  easily  construct 
an  equivalent  initial  value  problem  starting  with  a  plane  wave  initiated  at  some  finite 
(negative)  time.  (This  is  equally  true  for  a  wave  incident  from  the  left  and  below. 
This  is  the  case  extensively  treated  in  de  Hoop  [1958].)  For  9C  <  7t/2 — incidence  from 
the  right — this  is  not  so  straightforward.  For  any  finite  negative  time,  the  total  field 
would  have  to  include  a  plane  wave  terminated  on  the  reflector  plus  its  companion 
reflected  wave,  certainly  more  complicated  than  the  “ideal”  we  propose  to  study.  It 
will  be  seen  below  how  the  “noncausal”  nature  of  the  problem  posed  here  for  this 
latter  choice  of  9C  manifests  itself  in  the  analysis. 

In  the  mathematical  formulation,  one  begins  b)'  assuming  an  incident  wave 

V[(x,z,u)  =  exp{  —  icupi [xcos#c  +  2sin0c]},  p\  —  \jc\ ,  (0.93) 

with  incidence  angle 

9,  =  z  -  9C.  (0.94) 


35 


Bleistein  and  Cohen 


Cagniard  Method,  Complex  Time 


Note  t hat  in  this  problem,  we  take  z  positive  upward  and  we  take  9  to  be  the 
polar  angle  measured  from  the  positive  x  axis.  This  interchanges  the  role  of  sin#  and 
cos  6  from  the  standard  geophysical  convention.  However,  this  makes  comparison 
with  other  solutions  in  the  literature  easier,  since  this  is  Sommerfeld’s  notation. 

The  total  solution  is  written  as 


vT(x,z,u)  =  vj(x,z,uj)  +  vs(x,z,lj). 


(0.95) 


Here  the  scattered  field  v$(x,  z,uj)  will  include  the  reflected  field,  the  diffracted  field 
and  the  negative  of  the  incident  field  in  the  geometrical  shadow.  With  this  formulation, 
Vs(x,  z,u)  is  an  outward  propagating  wave  for  u  real.  Alternatively,  vs(x,  z,u)  must 
be  an  attenuating  wave  for  Imw  >  0,  or  equivalently,  Re  s  >  0,  which  is  where  the 
type  of  Fourier  transform  we  are  considering  here  is  initially  defined. 


The  solution  to  this  problem  by  the  Weiner-Hopf  technique  (Carrier,  Krook  and 
Pearson,  1983)  leads  to  the  integral  representation 


vs(x 


\  -  COS  6C)  f  eXp{*  [kx  +  ~k2\Z\  }  jf 

^  2ni  Jrt  [fc  +  up-,  ■  ■•'z9e\  :  r 


(0.96) 


The  contour  T*.  is  shown  in  Figure  0.6  for  the  case  Rea;  >  0  and  cos#c  >  0.  This  is 
the  case  which  is  inherently  noncausal. 


Recall  that  part  of  our  process  of  transformation  from  (k, u)  to  (p,s)  involved  a 
rotation  of  u  from  the  real  axis  to  the  imaginary  axis,  albeit  after  rescaling  k.  In 
Figure  0.6,  one  can  see  that  a  rotation  of  u  will  cause  the  pole  to  rotate  through  the 
path  of  integration.  This  is  not  acceptable,  unless  we  account  for  the  residue  at  the 
pole.  If  we  were  to  rescale  k ,  the  singularities  would  lie  on  the  axis;  for  Ima;  >  0  the 
image  contour  would  pass  between  them  from  upper  left  to  lower  right,  and  now  the 
same  difficulty  would  arise  from  further  rotation  of  the  contour  with  a  fixed  pole. 
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CO/C, 


-cocosGc/^ 


r. 


-CO  /c 


1 


Fig.  0.6.  The  contour,  T;.. 


When  cos Br  <  0,  the  pole  moves  to  the  first  quadrant.  IT  can  Ire  replace  by  the 
Re  A'  axis,  and  the  pole  rotates  away  from  the  contour  along  with  the  branch  point 
in  the  first  quadrant  as  argu?  increases.  Equivalently,  after  rescaling,  -k  —  u jp,  the 
contour  of  integration  would  move  away  from  all  singularities  with  increasing  a  just 
as  in  the  previous  two  examples.  Recall  that  this  is  the  case  that  corresponds  to  an 
equivalent  causal  problem  and  our  analysis  is  no  more  difficult  than  it  was  for  the 
standard  problems  treated  by  the  Cagniard-de  Hoop  method.  We  proceed  with  the 
analysis  of  that  more  interesting  cruse  for  which  there  is  no  equivalent  causal  problem, 
that  is,  for  cos 0r  >  0. 

Let  us  replace  T*.  in  Figure  0.6  by  the  contour  in  Figure  0.7  and  replace  the  integral 
over  this  contour  by  the  sum  of  the  residue  at  the  pole  and  an  integral  along  the  real 
axis.  That  is. 


vs(x,  =  vres(x,  +  v(x,  z,oj), 
i'res(x,z<u)  =  -  exp  {iojpi  [-xcos^c  +  |c|sin^c]}  ,  (0.97) 
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-cocosBc/^ 


-CO/C-I 


Fig.  0.7.  Replacement  contour  for  T^.. 


r(.r,  c.o.-)  =  - 


V^'FiU  -  cos Bc) 


o~, 


exp  |  i  kx  +  v/uAp'i  -  A'2|c|  i 

/  - !: _ _  1-Lflb 

J-oc  [A:  +  up i  cos#r]  \/ujp\  +  k 


Now  we  can  introduce  the  change  of  variables  (0.20),  k  —  u;/>sgn(.r),  carry  out 
the  rotation  of  lj  to  the  imaginary  axis,  and  set  u>  =  is  (as  previously  described)  .o 
obtain  the  standard  form  of  the  integral  representation  for  u(x,z,s)  —  v(x,z,is): 

/( 1  —  cos 6C )  f  exp{  —  sw{p.  |.r|,  c)} 


U(x,z,s)=  V 
Again,  in  this  conation. 


i  -  cos  acj  r  , 

2~iyJF\  Jr  sgn(,r 


[}>  sgn(.r)  +  Pi  cos 6C]  yjpsgn(x)  + 


--dp 


(0.98) 


Pi 


w 


(p,  |.r|,-)  =  p\.r\  +  \Jp'i  —  p2|~|,  p\  =  1/ci. 


(0.99) 


and  T  is  the  contour  of  Figure  0.2. 

This  analysis  could  have  started  from  an  integral  in  which  jj  was  real,  with  the 
contour  T*.  replaced  by  the  contour  shown  in  Figure  0.8.  The  final  result  would  not 
change.  The  rotation  of  contour  after  rescaling  k  and  rotation  of  u  or  the  rotation  of 
first  would  produce  exactly  the  same  interplay  between  the  contour  of  integration 
and  the  pole. 
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-CDCOS0C  /CA  CD/C-. 


— CO/C-, 


Fig.  0.8.  Contour  of  integration  for  Recu. 

We  have  now  made  the  main  point  of  this  example  concerning  the  Cagniard-de 
Hoop  method,  that  the  change  of  variables  plus  rotation  implicit  in  the  equation 
io  —  is ,  starting  from  u  (nearly)  real  and  ending  with  s  real,  may  introduce  special 
contributions  as  a  consequence  of  singularities  forced  to  pass  “through”  the  contour 
of  integration. 

It  is  interesting  to  continue  the  analysis  of  the  solution  via  the  Cagniard-de  Hoop 
method  and  we  proceed  to  do  so,  despite  having  already  accomplished  our  main 
objective  with  respect  to  this  problem. 

We  note  first  that,  vRES(x,  z,lj)  in  (0.97)  can  also  be  written  as 

{-  exp  (fcupi  [x  cos(n  -  9C)  -r  csin(7t  -  0C)]}  .  z  >  0, 

(0.100) 

-  exp  {iupi  [.tcos(7t  +  6C )  4-  zsin(7r  +  0C)]}  ,  c  <  0. 

The  first  line  here  has  the  form  of  the  reflected  wave,  while  the  second  line  has  the 
form  of  the  negative  of  the  incident  wave.  As  previously  noted,  we  are  alert  for  such 
terms  to  appear  in  the  solution  v$.  Ti  e  only  problem  is  that  these  two  expressions 
cover  the  entire  upper  and  lower  half-spaces,  respectively,  rather  than  the  reflection 
region  (labeled  I  in  Figure  0.9)  in  the  former  case  and  the  geometrical  shadow  (labeled 
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Fig.  0.9.  Geometrical  optics  regions  for  the  scattering  of  a  plane  wave  by  a  half-plane 
reflector.  Region  /:  geometrical  reflection.  Region  IT.  geometrical  shadow  for  the 
incident  wave.  Region  III :  remainder  of  the  plane,  v^es  extends  the  reflected  and 
negative  incidence  fields,  respectively,  to  Region  III. 

II  in  Figure  0.9)  in  the  latter  case.  Thus,  we  must  be  alert  to  further  corrections  to 
these  plane  waves  in  the  representation  of  u(x,z,s)  in  (0.98). 

We  now  turn  to  the  analysis  of  u(x,  z ,  s )  defined  by  (0.98).  Because  0  is  measured 
from  the  x-axis  in  this  example,  the  location  of  the  saddle  point  in  (0.34)  is  given  by 

/;«  =  pi  |  cos  0\  =  ■-°S^  =  — ,  p  =  Vx2  +  z2.  (0.101) 

Ci  C\p 

What  is  of  interest  to  us  here,  is  the  interplay  of  the  location  of  p3  and  the  pole  of 
the  integrand  in  (0.98).  We  note  first  that  the  pole  is  at  some  p  <  0  for  x  >  0  and 
therefore  is  outside  of  the  primary  p-domain.  Therefore,  we  need  only  consider  the 
case  x  <  0  and  the  pole  located  at 

pp  —  -sgn (x)pi  cos  6C  =  -sgn(x)cos^c/ci.  (0.102) 
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Fig.  0.10.  The  regions  of  incidence  and  reflection  described  in  terms  of  angular 
boundaries. 

In  this  case,  it  is  at  least  necessary  for  x  to  be  negative  for  the  pole  to  be  in  the 
primary  p-domain;  that  is, 

x  <  0  =s>  7r/2  <  6  <  37r/2,  | cos#|  =  —  cos#.  (0.103) 

Furthermore,  it  is  necessary  for  the  pole  to  be  to  the  left  of  the  saddle  point  for  it  to 
be  in  the  primary  p-domain;  that  is, 

Pp  <  P,  =>  cos  $c  <  —  cos  9  —  cos( 7r  —  9).  (0.104) 

When  we  rewrite  this  result  in  terms  of  angles, 

—  9C  <  n  -  6  <  0C  =s>  ir  —  6C  <  6  <  tt  +  9C,  (0.105) 

we  see  that  the  pole  is  in  the  primary  p-domain  for  (x,  z)  in  Region  III  of  Figure  0.9. 
See  Figure  0.10. 

For  this  angular  range,  the  image  of  the  pole  is  in  the  primary  r-domain  and 
the  deformation  of  T'  is  a  loop  integral  starting  at  the  pole  just  as  in  Figure  0.4. 
For  the  pole  in  the  r  plane  and  the  branch  point  at  r,  “well-separated,”  we  can 
consider  the  two  contributions  separately.  It  is  fairly  straightforward  to  verify  that 
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the  residue  at  this  pole  is  just  the  negative  of  the  wavefields  given  in  (0.100).  That 
is,  the  residue  at  this  pole  has  the  effect  of  eliminating  the  reflected  wave  and  the 
negative  of  the  incident  wave  in  the  regions  where  we  did  not  expect  them  to  appear 
in  the  first  place;  now  these  waves  are  restricted  to  the  region  of  reflection  and  to  the 
geometrical  shadow,  respectively. 


For  (z,  r)  near  the  boundaries  of  Region  III  the  pole  and  branch  point  are  nearby 
one  another.  If  we  wish  to  allow  them  to  coalesce,  we  cannot  consider  the  contribu¬ 
tions  separately.  In  the  limit  of  coalescence  of  the  singularities  we  no  longer  have  a 
simple  residue  plus  branch  cut  integral  and  we  must  content  ourselves  with  numerical 
integration  or  a  uniformly  valid  asymptotic  approximation  that  allows  coalescence. 


Let  us  suppose  now  that  the  pole  and  the  saddle  point  are  “well-separated,” 
whether  or  not  we  are  in  Region  III.  We  denote  by  uD(x,  y,  s)  the  contribution  from 
the  loop  integral  starting  at  the  saddle  point,  the  contour  S'  in  Figure  0.2.  This 
integral  is  obtained  in  the  standard  way,  as  described  for  the  previous  two  examples, 
from  the  representation  (0.98).  That  is, 


uD(x,z,s)  - 


-  cos 6C ) 


/, 


exp{— sr} 


dp 


'Iziyfcl  Js1  [psgn(z)  -j-  px  cos#c]  y'psgn(z)  +  p\  d T 


dr  .  (0.106) 


Since  w  in  (0.99)  is  the  same  as  (0.52),  the  analysis  of  the  first  example  can  be  used 
here,  as  well.  In  particular,  dp /dr  is  given  by  (0.61)  and  t,  is  given  by  (0.56).  After 
some  algebra,  we  find  in  analogy  with  (0.57)  that 


uD(x,z,s ) 


yj{  1  -  cos 0C)  too  e-sr 

niy/elp2  JPlp  ^r2  _  p2p2 


•Re 


T\Z\  —  l\X 


1/ 


r2  —  pfp2 


[p+(t)s gn(x)  +  Pl  cos6c]yJp+(r )sgn ( x )  -F  px 


|  dr.(0. 
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In  this  equation,  />+(r)  is  given  by  (0.62).  From  this  result,  the  time-domain  solution 
is  given  by 


uD{x,z,s) 


)/(l~cos  Bc)  H(t  —  pjc\ 
-K'sfclP1  ^Jt2  -  1/cfp2 


•Re 


t\z\  -  i\x\^T2  -  p2/c\ 


(0.108) 


b+(r)s§n(:r)  +  cosdc/ci]^/p+(r)sgn(r)  +  l/cx 
This  is  the  diffracted  wave.  It  is  a  cylindrical  wave  radiating  from  the  origin  at  time 
t  =  0,  which  is  the  arrival  time  of  the  incident  wave. 


We  have  seen  here  the  application  of  the  Cagniard-de  Hoop  method  to  a  “nontra- 
ditional”  problem.  The  transformation  from  (k,u)  to  (p,  s)  was  seen  to  have  a  new 
feature — the  interplay  of  the  rotation  of  contour  with  a  singularity  (a  pole)  of  the 
integrand.  Also,  for  this  problem,  the  body  wave  manifests  itself  through  a  pole  of 
the  integrand  and  the  saddle  point  of  w(p,x,  z)  gave  rise  to  the  edge-diffracted  wave. 
This  is  in  contrast  with  layered-media  problems  where  the  saddle  point  gives  rise  to 
the  body  wave. 


CONCLUSIONS 

We  have  proposed  an  approach  to  the  Cagniard  method  using  a  complex  time 
domain.  We  have  attempted  to  show  that  there  are  advantages  to  promoting  the 
time  domain  to  such  “equal  status”  with  the  slowness  domain.  We  avoid  the  usual 
folding  of  contour  in  the  slowness  domain  because  we  also  believe  that  there  are 
advantages  to  having  both  endpoints  of  integration  at  infinity.  Deformation  of  the 
original  contour  of  integration  onto  the  Cagniard  path  on  which  the  exponent  is  real 
is  equivalent  to  closing  down  the  image  contour  in  the  r-domain  on  the  real  axis. 
It  was  shown  in  this  process  that  a  particular  saddle  point  and  its  associated  paths 
of  steepest  ascent  play  a  crucial  role  in  the  analysis.  The  original  contour  and  the 
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contour  made  up  of  these  two  ascent  paths  define  a  primary  p-domain  whose  image, 
the  primary  r-domain,  is  central  to  the  analysis  of  deformation  of  contour.  Indeed, 
the  pair  of  steepest  ascent  paths  will  always  form  the  tails  of  the  usual  Cagniard 
paths  in  the  p-domain  and  will  make  up  those  paths  in  their  entirety  when  there 
are  no  other  singularities  of  the  integrand  in  the  primary  p-domain.  When  there  are 
singularities  in  the  primary  p-domain,  their  images  appear  in  the  primary  r-domain 
and  the  closing  down  of  the  image  contour  around  the  real  axis  in  that  domain  must 
take  account  of  these  singularities. 

We  have  also  addressed  the  problem  of  transforming  an  integral  solution  in  the 
( k,u )  domain  into  an  integral  solution  in  the  (p,  s)  domain.  It  was  shown  here  that 
careful  analysis  allows  us  to  address  noncausal  problems  as  well  as  the  causal  ones 
normally  treated  by  starting  with  one-sided  Laplace  transforms. 
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